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O ' Abstract 

. We provide an algebraic setting for cumulants and factorial moments through 

the classical umbral calculus. Main tools are the compositional inverse of the 
unity umbra, connected with the logarithmic power series, and a new umbra 
here introduced, the singleton umbra. Various formulae are given expressing 
. cumulants, factorial moments and central moments by umbral functions. 
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1 Introduction 

I ■ 
> 

I The purpose of this paper is mostly to show how the classical umbral calculus gives 

' a lithe algebraic setting in handling cumulants and factorial moments. The classical 

. umbral calculus consists of a symbolic technique dealing with sequences of numbers 

a„ indexed by nonnegative integers n = 0, 1, 2, 3, . . . , where the subscripts are treated 
I as if they were powers. This kind of device was extensively used since the nineteenth 

century although the mathematical community was sceptic of it, owing to its lack 



of foundation. To the best of our knowledge, the method was first proposed by Rev. 
John Blissard in a series of papers as from 1861 (cf. [5] for the full list of papers), 
nevertheless it is impossible to put the credit of the original idea down to him since 
the Blissard's calculus has its mathematical source in symbolic differentiation. In 
^ . the thirties, Bell [1] reviewed the whole subject in several papers, restoring the 

^ I purport of the Blissard's idea and in [2] he tried to give a rigorous foundation of the 

mystery at the ground of the umbral calculus but his attempt did not have a hold. 
Indeed, in the first modern textbook of combinatorics [12] Riordan often employed 
this symbolic method without giving any formal justification. It was first Gian- 
Carlo Rota to disclose the "umbral magic art" of shifting from o" to a„ bringing 
to the light the underlying linear functional (cf. [15]). This idea led Rota and his 
collaborators to conceive a beautiful theory (cf. [10] and [14]) which has originated 
a large variety of applications (see [4] for a list of papers updated to 2000). Some 
years later, Roman and Rota [13] gave rigorous form to the umbral tricks in the 
setting of Hopf algebra (see also [9]). But in 1994 Rota himself wrote (cf. [18]): "... 
Although the notation of Hopf algebra satisfied the most ardent advocate of spic-and-span 
rigor, the translation of "classical" umbral calculus into the newly found rigorous language 
made the method altogether unwieldy and unmanageable. Not only was the eerie feeling 
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of witchcraft lost in the translation, but, after such a translation, the use of calculus to 
simplify computation and sharpen our intuition was lost by the wayside..." Then, in the 
paper [18] The Classical Umbral Calculus (1994) Rota, together with Taylor, tries 
to restore the feeling meant by the founders of the umbral calculus keeping new 
notation both minimal and indispensable to avoid the misunderstanding of the past. 
In this new setting, the basic device is to represent an unital sequence of numbers by 
a symbol a, named umbra, i.e. to associate the sequence 1, ai, 02, . . . to the sequence 
1, a, a^, . . . of powers of a. through an operator E that looks like the expectation of 
random variables (r.v.'s). This new way of dealing with sequences of numbers has 
been applied to combinatorial and algebraic subjects (cf. [17], [23] and [8]), wavelet 
theory (cf. [19]) and difference equations (cf. [24]). Besides it has led to a nimble 
language for r.v.'s theory, as showed in [16] and [5]. 

The present work is inspired by this last point of view. As a matter of fact, 
an umbra looks as the framework of a random variable (r.v.) with no reference to 
any probability space, someway getting closer to statistical methods. However, the 
use of symbolic methods in statistics is not a novelty: for instance Stuart and Ord 
[22] resort to such a technique in handling moments about a point. In addition in 
the umbral calculus, questions as convergence of series are no matter, as showed 
hereafter dealing with cumulants. 

Among the sequences of numbers related to r.v.'s, cumulants play a central role 
characterizing all r.v.'s occurring in the classical stochastic processes. For instance, 
a r.v. having Poisson distribution of parameter x is the unique probability distri- 
bution for which all its cumulants are equal to x. It seems therefore that a r.v. 
is better described by its cumulants than by its moments. Moreover, due to their 
properties of additivity and invariance under translation, the cumulants are not nec- 
essarily connected with the moments of any probability distribution. We can define 
cumulants Kj of any sequence a„, n = 1, 2, 3, ... by 



in disregard of questions of whether any series converges. By this approach, many 
difficulties connected to the "problem of cumulants" smooth out, where with "prob- 
lem of cumulants" we refer to characterizations of sequences that are cumulants of 
some probability distributions. The simplest example is that the second cumulant 
of a probability distribution must always be nonnegative, and is zero only if all of 
the higher cumulants are zero. Cumulants are subject to no such constraints when 
they are analyzed by an algebraic point of view. What is more, in statistics they 
do not play any dual role compared to factorial moments. Whereas the algebraic 
setting here proposed comes to the light their close relationship through an umbral 
analogy with the well known complementary notions of compound and randomized 
Poisson r.v.'s (cf. [6]) 

Umbral notations are introduced in Section 2 by means of r.v.'s semantics. Our 
intention by this way is to make the reader comfortable with the umbral system of 
calculation without require any prior knowledge. We only skip some technical proofs 
of formal matters on which the reader is referred to citations. We also resume the 
theory of Bell umbrae, completely developed in [5], that not only gives the umbral 
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counterpart of the family of Poisson r.v.'s but it allows an umbral expression of the 
functional composition of exponential power scries. Section 3 is devoted to a new 
umbra, named singleton umbra, playing a dual role compared to the Bell umbra. 
Their relationship is encoded by the compositional inverse of the unity umbra. The 
singleton umbra is the keystone of the umbral presentation of cumulants and factorial 
moments. 

In the last two sections we give various umbral formulae for cumulants and 
factorial moments that parallel those known in statistics but simplify the proofs as 
well as the forms. This happens for instance for the equations expressing cumulants 
in terms of moments (and vice-versa) and also for their recursive formulas. Inversion 
theorems allowing to obtain an umbra from its cumulants or factorial moments are 
also stated. 

In 1929, Fisher [7] introduced the A;— statistics as new symmetric functions of the 
random sample. The aim of Fisher was to estimate the cumulants without using 
the moment estimators. He used only combinatorial methods. The A;— statistics are 
related to the power sum symmetric functions whose variables are the r.v.'s of the 
sample, but these expressions are very unhandy. We believe that the umbral calculus 
may seek to simplify the expression of the fc— statistics (as well as the /i— statistics 
for the central moments) taking into account its combinatorial nature. 

2 Umbrae and random variables 

In the following, we resume terminology, notations and some basic definitions of 
the classical umbral calculus, as it has been introduced by Rota and Taylor in [18] 
and further developed in [5]. Fundamental is the idea of associating a sequence of 
numbers 1, 0,2, a^, . . . to an indeterminate a which is said to represent the sequence. 
This device is familiar in probability when represents the i— th moment of a 
r.v. X. In this case, the sequence 1, ai, 02, . . . results from applying the expectation 
operator E to the sequence 1, X, X^, . . . consisting of powers of the r.v. X. 
More formal, an umbral calculus consists of the following data: 

a) a set A = {a, f3, . . .}, called the alphabet, whose elements are named umbrae; 

b) a commutative integral domain R whose quotient field is of characteristic zero; 

c) a linear functional E, called evaluation, defined on the polynomial ring R[A] and 

taking values in R such that 

i) E[l] = 1; 

ii) E[q^(3^ ■ ■ •7'^'] = E[d'-\E[(3^ ■ ■ ■ E[^^] for any set of distinct umbrae in A 

and for i,j,...,k nonnegative integers (uncorrelation property); 

d) an element e ^ A, called augmentation [13], such that £^[e"'] = (5o,n) for any 

nonnegative integer n, where 

I 1 if i = 7 
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an element u E A, called unity umbra [5], such that Elu^] = 1, for any nonnega- 
tive integer n. 



A sequence ao = 1, ai, 02, ■ ■ ■ in -R is umbrally represented by an umbra a when 

E[a'] = ai, for i = 0,1,2,... . 

The elements Oj are called moments of the umbra a on the analogy of r.v.'s theory. 
The umbra e can be view as the r.v. which takes the value with probability 1 and 
the umbra u as the r.v. which takes the value 1 with probability 1. Note that the 
uncorrelation property among umbrae parallels the analogue one for r.v.'s as well 
as it is E[a"''^'^] ^ E[a'^]E[a'']. Remark as this setting gets out of the well-known 
"moment problem" for r.v.'s. 

Example 2.1 Bell umbra. 

The Bell umbra (3 is the umbra such that 

E[{p)^] = l n = 0,1,2,... 

where (/3)o = 1 and (/9)„ = — 1) • • • (/3 — n + 1) is the lower factorial. It results 
= Bn where Bn is the n— th Bell number (cf. [5]), i.e. the number of the 
partitions of a finite nonempty set with n elements or the n— th coefficient in the 
Taylor series expansion of the function exp(e* — 1). So /? is the umbral counterpart 
of the Poisson r.v. with parameter 1. 

We call factorial moments of an umbra a the elements 

_ f 1, n = 

"(")-\ E[{aU n>0 

where (a)„ = a{a — 1) • • • (a — n + 1) is the lower factorial. So the definition of P 
in example 2.1 could be reformulated as follows: the Bell scalar umbra is the umbra 
whose factorial moments are = 1 for any nonnegative integer n. 



2.1 Similar umbrae and dot-product 

The notion of similarity among umbrae comes in handy in order to manipulate 
sequences such 

E " CLiGn-i, ne N (1) 

as moments of umbrae. The sequence (1) cannot be represented by using only the 
umbra a with moments oq = 1, 04, 0,2^ • • • • Indeed, being a correlated to itself, the 
product aiUn-i cannot be written as E[a'''a'^~^]. So we need two distinct umbrae 
having the same sequence of moments, as it happens for similar r.v.'s. Therefore, if 
we choose an umbra a' uncorrelated with a but with the same sequence of moments, 
it is 

E " aiO-n-i = E 
i=o V W 




a' {a') 



l\n- 



E[{a + a 



(2) 
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Then the sequence (1) represents the moments of the umbra (a + a'). A way to 
formahzc this matter is to define two equivalence relations among umbrae. 
Two umbrae a and 7 are umhrally equivalent when 



in symbols = 7. We note that equality implies similarity which implies umbral 
equivalence. The converses are false. Then, we shall denote by the symbol n.a the 
dot-product of n and a, an auxiliary umbra (cf. [18]) similar to the sum a' + a" + 
. . . + a'" where j are a set of n distinct umbrae each similar to the 

umbra a. So the sequence in (2) is umbrally represented by the umbra 2. a. We 
assume that O.a is an umbra similar to the augmentation e. 

We shall hereafter consider the dot product of n and a as an umbra if we saturate 
the alphabet A with sufficiently many umbrae similar to any expression whatever. 
For a formal definition of a saturated umbral calculus see [18]. It can be shown that 
saturated umbral calculi exist and that every umbral calculus can be embedded in 
a saturated umbral calculus. 

The following statements are easily to be proved: 

Proposition 2.2 (i) If n.a = n.p for some integer n ^ then a = /?; 

(ii) if c £ R then n.{ca) = c{n.a) for any nonnegative integer n; 

(Hi) n.{m.a) = (nm).a = m.(n.a) for any two nonnegative integers n,m; 

(iv) (n + m).a = n.a + m.a' for any two nonnegative integers n,m and any two 

distinct umbrae a = a'; 

(v) {n.a + n.l3) = n.[a + (3) for any nonnegative integer n and any two distinct 

umbrae a and (3. 

Two umbrae a and 7 are said to be inverse to each other when a + 7 = e. We 
denote the inverse of the umbra a by — l.a', with a = a' . Recall that, in dealing 
with a saturated umbral calculus, the inverse of an umbra is not unique, but any 
two inverse umbrae of the same umbra are similar. 

Example 2.3 Uniform umbra. 

The Bernoulli umbra (cf. [18]) represents the sequence of Bernoulli numbers S„ 
such that 



The inverse of the Bernoulli umbra is the umbral counterpart of the uniform r.v. 
over the interval [0,1] (cf. [23]). 



E[a\=E[^], 



in symbols a c± p. They are similar when 



a"~7", n = 0,1,2,... 
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2.2 Generating functions 

The formal power series in 



is the generating function (g.f.) of the umbra a, and it is denoted by e"*. The notion 
of umbrally equivalence and similarity can be extended coefficientwise to formal 
power series i?[yl][[i]] (see [24] for a formal construction). So it results 

a = /3 ^ e"* ~ e^*. 

Moreover, any exponential formal power series^ in R[[t]] 

ffi 

n>l 

can be umbrally represented by a formal power series (3) in In fact, if the 

sequence l,oi,a2, . . . is umbrally represented by a then 

f{t) = E[e'''] i.e. /(t):^e"*, 

assuming that we naturally extend E to be linear. We will say that f{t) is umbrally 
represented by a. Note that, from now on, when there is no mistaking, we will just 
say that f{t) is the g.f. of a. For example the g.f. of the augmentation umbra e is 1 
as well as the g.f. of the unity umbra u is e^. 

Getting back to a r.v. X, recall that when E[exp(tX)] is a convergent function 
f(t), it admits an exponential expansion in terms of the moments which are com- 
pletely determined by the related distribution function (and vice- versa). In this case 
the moment generating function (m.g.f.) encodes all the information of X and the 
notion of similarity among r.v.'s corresponds to that of umbrae. 

The first advantage of the umbral notation introduced for g.f.'s is the represen- 
tation of operations among g.f.'s with operations among umbrae. For example the 
multiplication among exponential g.f.'s is umbrally represented by a summation of 
the corresponding umbrae: 

5(t)/(t) :^ e("+^)* with /(t) e-*, 3(t) =^ e^*. (4) 

Via (4), the g.f. of n.a is f{t)^. If a is an umbra with g.f. f{t), the inverse um- 
bra —l.a' has g.f. [/{t)]"^. The summation among exponential g.f.'s is umbrally 
represented by a disjoint sum of umbrae. The disjoint sum (respectively disjoint 
difference) of a and 7 is the umbra 77 (respectively t) with moments 

n ] u, n = ( ^. , n \ u, n = 0\ 

^ - ^ n , n ^ A respectively t ~ < „ „ ^ n ' 
1 a" -I- 7", n > V 1 a" - 7", n > / 

in symbols 77 = a-i-7 (respectively l = a— 7). By the definition, it follows 

f{t) ± [g{t) - 1] ~ e("±^)*. 



^Observe that with this approach we disregard of questions of whether any series converges. 
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Example 2.4 Unbiased estimators. 

Suppose to make the disjoint sum of n times the umbra a. We will denote this umbra 
by Its g.f. is 1 + n[f{t) — 1]. The umbra +„« has the following probabilis- 

tic counterpart. Let {Xj}"^^ be a random sample of independent and identically 
distributed (i.i.d.) r.v's. As it is well-known the power sum symmetric functions 

n 
i=l 

gives the unbiased estimators Sr/n of the moments of Xi. But E[{+nC(y] = nur, 
hence the umbral corresponding of the power sum symmetric functions sequence Sr 
is the umbra -i-nO;- 

2.3 Auxiliciry umbrae 

In the following, suppose a an umbra with g.f. f{t) and 7 an umbra with g.f. g{t). 
The introduction of the g.f. device leads to the definition of new auxiliary umbrae 
useful for the development of the system of calculation. For this purpose, we should 
replace R with whatever polynomial ring having coefficients in R and a number 
of indeterminates according to necessity. In this paper, we deal with R[x,y]. This 
allows to define the dot-product of x and a via g.f., i.e. x.a is the auxiliary umbra 
having generating function 

The Proposition 2.2 still holds replacing n with x and m with y. Then, an umbra is 
said to be scalar if the moments are elements of R while it is said to be polynomial 
if the moments are polynomials. 

Example 2.5 Bell polynomial umbra. 

The Bell polynomial umbra (j) is the umbra having factorial moments equal to x'^ (cf. 
[5]). This umbra has g.f. exp[x(e* — 1)] so that <p = x.p, where P is the Bell umbra. 
It turns out that the Bell polynomial umbra x.p is the umbral counterpart of the 
Poisson r.v. with parameter x. 

Example 2.6 Moments about a point. 

The moments E[{X — a)"] about a point a G R of a r.v. X are easily represented 
by umbrae through the following definition: the umbra a" having moments about a 
point a € R is defined as 

a" = q; — a.u. (5) 

If a, 6 G i? and b — a = c, then 

a"' = a — (b + c).u = a'^ + c.u, 

is the umbral version of the equations giving the moments about a in terms of the 
moments about b (cf. [22] for another symbolic expression). 

The dot-product 7.0 of two umbrae is the auxiliary umbra having g.f. 

e(7.«)t ^ [/(^)]7 ^ e7iog/(t) ~ ^[iog/(t)] . 
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The moments of the dot-product 7.0; are (cf. [5]) 

n 

E[{'^.ay-]=^g{i)Bn,i{ai,a2,...,an-i+i) n = 0,1,2,... (6) 

1=0 

where (7(j) are the factorial moments of the umbra 7, Bn^i are the (partial) Bell 
exponential polynomials (cf. [12]) and are the moments of the umbra a. Observe 
that £^[7.0;] = giai= -£[7] E[a] The following properties hold (cf. [5]): 

Proposition 2.7 a) ifrj.a = rj.^ then a = j; 

b) if c E R then rj.{ca) = 0(77.0;) for any two distinct umbrae a and rj; 

c) if J = 7' then {a + r]).j = a.7 + 77.7'; 

d) r].{'j.a) = {ri.^).a. 

Observe that from property b) it follows 

a.x = a.{xu) = x{a.u) = xa. (7) 

Remark 1 The auxiliary umbra 7.0 is the umbral version of a random sum. Indeed 

the m.g.f. ^[log f{t)] corresponds to the r.v. Sn = Xi + X2-\ \-Xn where iV is a 

discrete r.v. having m.g.f. g{t) and Xi are i.i.d. r.v.'s having m.g.f. f{t). The right- 
distributive property of the dot-product 7.0; runs in parallel with the probability 
theory because the random sum Sn+m is similar to Sn + Sm, where N and M are 
independent discrete r.v.'s. The left-distributive property of the dot-product 7.0 
does not hold as well as it happens in the r.v.'s theory. In fact, let Z = X + Y he 
a r.v. with X and Y independent r.v.'s. As it is easy to verify, the random sum 
Sn = Zi + Z2 + • ■ ■ + Zn, with Zi i.i.d. r.v.'s similar to Z, is not similar to the r.v. 
Sn+S% with S§ =Xi+X2 + --- + Xn, and Xi i.i.d. r.v.'s similar to X and with 
S^ = Yi+Y2-\ \-Yn and Yi i.i.d. r.v.'s similar to Y. 

Example 2.8 Randomized Poisson r.v. 

Let us consider the Bell polynomial umbra x.p. If in the place of x we put a generic 
umbra a, we get the auxiliary umbra a./? whose factorial moments are 

(a./3)„~a" ra = 0,1,2,... 

and moments given by the exponential umbral polynomials (cf. [5]) 

n 

(a./?)'' ~ $n(a) - n = 0,1,2,.... (8) 

i=0 

Its g.f. is /[e* — 1]. The umbra a. (3 represents a random sum of independent Poisson 
r.v.'s with parameter 1 indexed by an integer r.v. Y, i.e. a randomized Poisson r.v. 
with parameter Y. 
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As suggested in [14] , there is a connection between compound Poisson processes and 
polynomial sequence of binomial type, i.e. sequence {pn{x)} of polynomials with 
degree n satisfying the identities 

Pn{x + y) =J2^^.^ Pi{x)Pn-i{y) 

for any n (cf. for instance [10]). Two different approaches can be found in [3] 
and in [21]. A natural device to make clear this connection is the a— partition 
umbra I3.a, introduced in [5]. Its g.f. is exp[/(t) — 1] and it suggests to interpret 
a partition umbra as a compound Poisson r.v. with parameter 1. As well-known, 
a compound Poisson r.v. with parameter 1 is introduced as a random sum Sn = 
X\ + X2 + • • • + Xjv where has a Poisson distribution with parameter 1. The 
umbra (3.a fits perfectly this probabilistic notion taking into consideration that the 
Bell scalar umbra (3 plays the role of a Poisson r.v. with parameter 1. What's 
more, since the Poisson r.v. with parameter x is umbrally represented by the Bell 
polynomial umbra a compound Poisson r.v. with parameter x is represented by 
the polynomial a— partition umbra x.ip = x.(3.a with g.f. cxp[,T(/(t) — 1)]. The name 
"partition umbra" has a probabilistic ground. Indeed the parameter of a Poisson r.v. 
is usually denoted hy x = Xt, with t representing a time interval, so that when this 
interval is partitioned into non-overlapping ones, their contributions are stochastic 
independent and add to Sn- This last circumstance is umbrally expressed by the 
relation 

{x + y).(3.a = x.(3.a + y.(3.a (9) 

giving the binomial property for the polynomial sequence represented by x./3.a. In 
terms of g.f.'s, the formula (9) means that 

h^x+y{t) = K{t)hy{t) (10) 

where hx{t) is the g.f. of x.l3.a. Viceversa every g.f. hx{t^ satisfying the equality 
(10) is the g.f. of a polynomial a— partition umbra. The a— partition umbra repre- 
sents the sequence of partition polynomials = y^(ai, 02, . . . , a„) (or complete Bell 
exponential polynomials [12]), i.e. 

n 

E[{p.a)"-] = ^Bn,i{ai,a2,...,an-i+i) = 1^1(01, 02, ctn), (H) 

i=0 

where are the moments of the umbra a. Moreover every a— partition umbra 
satisfies the relation 

(/3.a)" ~ «'(/?.« + a')"~^ a^a', n = 0,1, 2,... (12) 

and conversely (see [5] for the proof). The previous property will allow an useful 
umbral characterization of the cumulant umbra (see corollary 4.12 in section 4.) The 
umbra p.a plays a central role also in the umbral representation of the composition of 
exponential g.f.'s. Indeed, the composition umbra of a and 7 is the umbra r = ^.(3. a. 
The umbra r has g.f. g[f{t) — 1] and moments 

n 

-^M = ^9iBn,i{(^i,a2, ■ ■ ■ ,an-i+i) (13) 

i=0 
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with gi and Oj moments of the umbra 7 and a respectively. We denote by a^"^^ the 
compositional inverse of a, i.e. the umbra having g.f. f~^{t) such that — = 
f[f~^{t) — 1] = 1 + t. For an intrinsic umbral expression of the compositional inverse 
umbra see [5], where it is also stated an umbral version of the Lagrange inversion 
formula. 

Example 2.9 Randomized compound Poisson r.v. 

As already underlined in example 2.8, the umbra 7./? represents a randomized Pois- 
son r.v. Hence it is natural to look at the composition umbra as a compound 

randomized Poisson r.v., i.e. a random sum indexed by a randomized Poisson r.v. 
Moreover, being (7./?). a = 7.(/3.a) (cf. statement d) of Proposition 2.7), the previ- 
ous relation allows to see this r.v. from another side: the umbra j.{p.a) generalizes 
the concept of a random sum of i.i.d. compound Poisson r.v. with parameter 1 in- 
dexed by an integer r.v. X, i.e. a randomized compound Poisson r.v. with random 
parameter X. 

At the end, the symbol a " denotes an auxiliary umbra similar to the product 
a' a" ■ ■ ■ a'" where a' , a" , . . . , a'" are a set of n distinct umbrae each similar to the 
umbra a. We assume that a'^ is an umbra similar to the unity umbra u. The mo- 
ments of a " are: 

£[(«-)'=] = £;[(«'=)•"] = a^, fe = 0,1,2,..., (14) 

i.e. the n— th power of the moments of the umbra a. Thanks to this notation in [5], 
the umbral expression of the Bell exponential polynomials was given as follows: 

— - -...(;)-.. » 

whenever oi 7^ and where a is the umbra with moments 
Example 2.10 The central umbra. 

We call central umbra the umbra a'^^ having moments about ai = E[a]. From (5), the 
classical relation between moments and central moments of a r.v. has the following 
umbral expression: 



(a«i)'* - E ( ? J i-lT~''{a')''a-^''-''\ a = a' n= 1,2,. 




k=0 



being E[{ai.u)''] = = E[a-^] from (14). 



3 The singleton umbra 

The singleton umbra plays a dual role compared to the Bell umbra, even if it has 
not a probabilistic counterpart. Besides, the singleton umbra turns out to be an 
effective symbolic tool in order to umbrally represent some well-known r.v.'s as well 
as cumulants and factorial moments. 
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Definition 3.1 (The singleton umbra) An umbra % is said to be a singleton um- 
bra if 

X"^^i,n n=l,2,---. 
The g.f. of the singleton umbra % is 1 + i- 
Example 3.2 Gamma r.v. 

The m.g.f. of a Gamma r.v. with parameters a and c is 



M{t) 



(1 - ct)° 



This g.f. is umbrally represented by the inverse of — c(a.x) (see (ii) of Proposition 
2.2 replacing nhy a E R). 

Table 1 lights up the duality between the singleton umbra x and the Bell umbra /?. 



Umbra 


Generating 


function 




X 


(i+t) = i+E;r=i 




x-X 


{l+tr = l + En=l 






a.x 


/[log(l + t)] 






l + log[/(t)] 






e^^-i = 1 + E^=i 






x.p 


^x{e*-l) ^ 1 + ^^^^ 


Yrk=iS{n,k)x^ 


t" 


a. (5 


/(e*- 


-1) 






e/(*) 


-1 





Table 1. 

The connection between the singleton umbra x the Bell umbra /3 is made 
clear in the following proposition. 

Proposition 3.3 Let x the singleton umbra,, fi the Bell umbra and u'^~^^ the 
compositional inverse of the unity umbra u. It results 



X = u<-^>.p = p.u<-^>, 
13.x = u = x-P- 



(17) 
(18) 



Proof. The g.f. of ti^ .f5.u is 1 + t, being u'" and u compositional inverses. 
So equivalence (17) follows by property a) of proposition 2.7 being 



u 



<-i> 



^X^X-'^ 



Equivalence (18) follows via g.f.'s in Table 1. ■ 

Distributive properties of the singleton umbra respect to the sum and the disjoint 
sum of umbrae are given in the following. 



Proposition 3.4 It results 



X-(a + 7) = X-a+X-1 
(a-i-7).X = a.x+7-X- 



(19) 
(20) 



11 



Proof. Let f{t) be the g.f. of a and g{t) the g.f. of 7. Equivalence (19) follows 

observing that the g.f. of x-(« + /9) is 1 + log[f{t)g{t)] = l + Iog[/(t)] +log[5((t)], i.e. 
the g.f. of x-Q^+X-/5- Equivalence (20) follows observing that the g.f. of {a+P)-X is 
/[log(l + t)] + 5[log(l + t)] - 1, i.e. the g.f. of a.x+P-X- ■ 

The notion of mixture of r.v.'s has an umbral counterpart in the disjoint sum +. 
Indeed let {ai}f^i be n umbrae and G -R be n weights such that 

n 
i=l 

The mixture umbra 7 of {aj}^^^ is the following weighted disjoint sum of {ai}f^i 
7 = x-Pi-P-ai+x-P2-P-a2+ ■ ■ ■ +X-Pn-(3-an (21) 

where /3 is the Bell umbra and % is the singleton umbra. Prom (19) equivalence (21) 
can be rewritten as 

7 = X-(Pi-/9-ai + P2-f3.a2 + • . . + p„./3.a„). 

Since the g.f. of J2'^=iPi-P-ai is expQ2i=iPi[fiit) - 1]), where fi{t) is the g.f. of a^, 
from Table 1 it follows that the g.f. of 7 is Y^i=iPifi{t). 

Example 3.5 Bernoulli umbral r.v. 

Let us consider the Bernoulli r.v. X of parameter p. Its m.g.f. is g{t) = q + pe* with 
q = 1 — p. The Bernoulli umbral r.v. is the mixture of the umbra e and the unity 
umbra u : 

C = X-Q-P-^+X-P-P-U- 

Recalling that x-Q-P-^ = £ it is 

^ = X-P-P- 

Indeed it is 

E[e^^] = 1 + logfe^'^^'-i)] =q + peK 
Example 3.6 Binomial umbral r.v. 

As it is well-known a binomial r.v. Y with parameters n G N, p G [Oj 1]) is the sum 
of n i.i.d. Bernoulli r.v.'s having parameter p. Then the binomial umbral r.v. is 

n.^ = n.x-P-P- 

The parallelism is evident if we recall that the m.g.f. of the binomial r.v. Y is 

fit) = iq + pe'r. 

4 The cumulant umbra 

For a r.v. having moments ai, 02, . . . , a„ and cumulants ni, K2, ■ ■ ■ , Hn it is 

On = ^ CttKtt and = X] ^■^"'■^ (^^) 

TT TT 
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the sums here are taken over the partitions tt = [j^^,j2^^, ■ ■ ■ of the integer n, 

and 



n 

C-jr 



0'i!)™i(j2!)"'^ • • • (Jfc0™'= milma! ■■■mk\ 

d-iT = CT,{-iy''~^{u.^ - ly. and i^^ = mi + m2 H h 

a,r = and = JJ^ Kj. 

In this section we show how the umbral calculus simplifies the above expressions, as 
well as the recursive formulae which give moments in terms of cumulants. 
Let a be an umbra with g.f. f{t). 

Definition 4.1 The cumulant of an umbra a is the umbra Ka defined by 

Ka = X-a 

where x is the singleton umbra. 

Definition 4.1 gives the umbral version of the second equality in (22). Moreover the 
first moment of the cumulant umbra Ha is ai, i.e. the first moment of the umbra a, 
being E[k^] = E[x] E[a] = E[a] = ai. 

Example 4.2 Cumulant of the umbra e. 

Since e = %.£, the umbra e is the cumulant umbra of itself, i.e. = e. 
Example 4.3 Cumulant of the umbra u. 

Since x = X-^^j the umbra x is the cumulant umbra of the umbra -u, i.e. Ku = X- 
Example 4.4 Cumulant of the Bell umbra. 

Since u = x-P (see (18)), the umbra u is the cumulant umbra of the Bell umbra P, 
i.e. Kfj = u. Prom example 2.1, the Poisson r.v. of parameter 1 has cumulants equal 
to 1. 

Proposition 4.5 The cumulant umbra Kq, has g.f. 

k{t) = l + \og[f{t)]. (23) 
Proof. See Table 1. ■ 
Example 4.6 Cumulant of the singleton umbra. 

Since 1 + log(l + t) is the g.f. of the umbra this umbra is the cumulant 

umbra of the umbra x, i.e. k,^ = 

Example 4.7 Cumulant of the Bernoulli umbral r.v. 

From example 3.5, the cumulant umbra of the Bernoulli umbral r.v. is X-iX-P-P) = 
u<-i>.p./3. 
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Example 4.8 Cumulant of the Binomial umbral r.v. 

Prom example 3.6, the cumulant umbra of the Binomial umbral r.v. is x.(n.x.p./3), 
i.e. 

t I ol < II oil < I aiii \ 

where x') x"; • • • > x'" S-^'S a set of n distinct umbrae each similar to the singleton 
umbra x as well as . . . ,/?'" are a set of n distinct umbrae each similar to the 

Bell umbra /3. Prom (19) and recalling examples 2.4 and 4.6, it results 

X-n.X-P-P = X-X'-P-P'+X-X"-P-P"+ ■ ■ ■ +X-x!"-P-0" = +nU<~^>.p.p. 
This parallels the analogous result in probability theory. 

From (23), the moments of the cumulant umbra 

{ka)n = EK] = 

that is equivalent to the definition of the n— th cumulant of a r.v. X having m.g.f. 
fit). 

To state the explicit version of the second equality in (22) 

n 

kn = J2{-^r\i - l)!S„,i(ai, a2, . . • ,a„_i+i) (24) 

i=l 

giving cumulants in terms of moments, usually requires laborious computations (cf. 
for example [11]). The umbral definition of cumulants allows a simple proof of (24). 
Indeed, being x = the cumulant umbra of a is the umbral composition of 

and a : 

and then its moments are given by (13). Equality (24) follows recalling that the 
moments of u'^"^^ are the coefficient of the exponential expansion 

oo ,i 

i + iog{i + t) = i + Y,{-iy-\i-i)\-. 

1=1 

Similarly, the three main algebraic properties of cumulants can be easily recovered 
from next theorem. 

Theorem 4.9 It is 

a) (the additivity property) 

X.{a + j) = x-a+x-l, (25) 

i.e. the cumulant umbra of a sum of two umbrae is equal to the disjoint sum 
of the two corresponding cumulant umbrae; 

b) (the semi-invariance under traslation property) for any c E R 

X.{a + c.u) = x-a+x-c; 
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c) (the homogeneity property) for any c E R 

X.(ca) = c(x.a). 

Proof. Property a) follows from (19). Property b) follows from (25), setting 7 = 
c.u for any c G -R. At the end, property c) follows from 6^ of proposition 2.7. ■ 

Example 4.10 Cumulant of the central umbra. 

The sequence of cumulants related to the central umbra a"^ is the same of a except- 
ing the first equal to 0. Indeed, by the additivity property of the cumulant umbra 
it is 

X.{a - ai.u) = x-a-x-ai- 

The results follows from (7). 

The umbral version of the first equality in (22) is given in the following theorem. 
Theorem 4.11 (Inversion theorem) Let Kq, be the cumulant umbra of a, then 

a = p.Ka 

where /? is the Bell umbra. 
Proof. It is 

p.Ka = /3.u'^~^^ ./3.a = X-/3-Q: = u.a = a. 



The inversion theorem allows to calculate the moments of the umbra a according to 
its cumulants. Recalling (11) it is 

an = Yn[{ka)l, {ka)2, ■ ■ {ka)n] (26) 

with a„ the n— th moment of the umbra a and {ka)n the n— th moment of the umbra 
Ka- Equation (26) is the explicit version of the first equality in (22). 

Remark 2 The complete Bell polynomials in (11) are a polynomial sequence of 
binomial type. Since from the inversion theorem any umbra a could be seen as 
the partition umbra of its cumulant Ka, it is possible to prove a more general result: 
every polynomial sequence of binomial type is completely determined by its sequence 
of formal cumulants. Indeed, in [5] it is proved that any polynomial sequence of 
binomial type represents the moments of a polynomial umbra x.a and viceversa. So 
from the inversion theorem any polynomial sequence of binomial type represents the 
moments of a polynomial umbra x.^.Ha- 

The next corollary follows from (12) and from the inversion theorem. 
Corollary 4.12 // the cumulant umbra of a, then 

a" ~Kc«(«a + «)"-' (27) 

for any nonnegative integer n. 
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Equivalences (27) were assumed by Shen and Rota in [16] as definition of the cumu- 
lant umbra. In terms of moments, equivalences (27) give 




that is largely used in statistic framework [20]. 
Example 4.13 Levy process. 

Let {Xt, t > 0) he a real- value Levy process, i.e. a process starting from and with 
stationary and independent increments. According to the Levy-Khintchine formula 
(cf. [6]), if we assume that Xt has a convergent m.g.f. in some neighbourhood of 0, 
it is 

E[e^^*] = 6*^=^ (28) 

where k{d) is the cumulant g.f. of X^. The inversion theorem gives the umbral 
version of equation (28): 

t.a = t.p.Ka- 
4.1 Cumulants of the Poisson r.v.'s 

From example 2.9, the umbra j.p.a corresponds to a compound randomized Poisson 
r.v., i.e. a random sum Sn = Xi + • • • + X^v with N a randomized Poisson r.v. of 
parameter the r.v. Y. In particular a corresponds to X and 7 corresponds to Y. 
Since x-{l'-P-<^) = the cumulant umbra of the composition of a and 7 is the 

composition of a and k^. Then from (13), the cumulants of a compound randomized 
Poisson r.v. are given by 

n 

^kiBn,i{a-i,a2,...,an-i+i) (29) 

i=l 

where a, are the moments of the r.v. X and ki are the cumulants of the r.v. Y. Now 
set 7 = x.u in 7./?. a. This means to consider a r.v. Y such that P{Y = x) = 1. Then, 
the random sum Sn becomes a compound Poisson r.v. of parameter x corresponding 
to the polynomial a— partition umbra x.p.a, with a the umbral counterpart of X. 
and cumulants 

n 

^kiBn,i{ai,a2,...,an-i+l) ~ XOn. (30) 

i=l 

Indeed (30) follows from (29) since the moments fcj of x-x are equal to 0, except 
the first equal to x. If x = 1 the cumulant of a— partition umbra is a so that the 
moments of X are the cumulants of the corresponding compound Poisson r.v. Now, 
in x.(3.a take a = u. From (30), the cumulants of the Bell polynomial umbra x.p 
are equals to x as well as for the Poisson r.v. of parameter x. 

At the end, in j.p.a set a = u. The cumulant umbra of 7./? is Kj.(3 with the 
cumulant umbra of 7. Its probabilistic counterpart is a randomized Poisson r.v. of 
parameter the r.v. Y, corresponding to the umbra 7. From (8) the cumulants of a 
randomized Poisson r.v. of parameter the r.v. Y are the moments of k^./3, i.e. 

n 

^S{n,i)ki 

1=0 
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with ki the cumulants of the r.v. Y. 



5 The factorial umbra 

The factorial moments of a r.v. do not play a very prominent role in statistics, but 
they provide very concise formulae for the moments of some discrete distributions, 
like the binomial one. 

Let a be an umbra with g.f. f{t). 

Definition 5.1 An umbra (pa is said to be an a— factorial umbra if 



where x is the singleton umbra. 
Example 5.2 e— factorial umbra. 

Since e = £.%, the e— factorial umbra is similar to the umbra e, i.e. ip^ = e. 
Example 5.3 u— factorial umbra. 

Since x = '^■X-, t^ie n— factorial umbra is similar to the umbra Xi i-^- 'Pu = X- 
Example 5.4 f3— factorial umbra. 

Since u = (3.x from (18), the /3— factorial umbra is similar to the unity umbra u, i.e. 
(/?^ = u. Prom example 2.1, the Poisson r.v. of parameter 1 has factorial moments 
equal to 1. 

Example 5.5 x~/<ictona/ umbra. 

Prom example 4.6, it is x-X = u'^~^^ ■ The %— factorial umbra turns out to be u"^"^^ , 
i.e. ip-^ = u"^"^^. 

Proposition 5.6 The a— factorial umbra has g.f. 



The a— factorial umbra has moments equal to the factorial moments of the umbra 
a, as the following proposition shows. 

Proposition 5.7 Let ipa be an a— factorial umbra. Then 



fa = a-X 



9{t) = f[\og{l+t)]. 



(31) 



Proof. See Table 1. 



iplc^{a)n,n = 0,1,2,.... 



Proof. By equation (6) and definition 5.1 it is 



n 




(32) 
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Since the umbra x moments equal to and x 



.k 



1, then 



Bn,k{h,\-, (^1,2, ■ ■ ■ , h,n-k+l) — 

Hence the equation (32) becomes E[{ipa)^] = (a) 



0, 
1, 



if n > k 
if n = k. 



(33) 



Example 5.8 Factorial umbra of the central umbra. 
Prom property c) of Proposition 2.7, it is 

= {a- ai.u).x = a.x - ai-x' = (fa - 



with x' = X- Then the factorial umbra of the central umbra a""^ is the difference 
between the factorial umbra of a and the factorial umbra of the umbra having 
moments equal to oi. By (31) its g.f. results /[log(l + t)]{l — t)"'^. 

Example 5.9 Factorial moments of the binomial r.v. 

Since the factorial moments characterize the binomial r.v., we show how to evaluate 
them by umbral methods. As showed in example 3.6, the umbral counterpart of the 
binomial r.v. is n.x-P-P- Due to (18) and (7) the corresponding factorial umbra is 
n-{x-P-P)-X = n-X-P = P{n.x)- Its g.f. is 



and so the factorial moments are {n)jp' . If n = 1, the factorial umbra is and 
from example 3.5 the first factorial moment of the Bernoulli r.v. is equal to p while 
the others are equal to 0. 

Example 5.10 Factorial umbra of the cumulant umbra. 

If Ka is the cumulant umbra of a, then Ka-X is the factorial cumulant umbra of a, 



The following theorem allows to obtain the umbra a from its factorial umbra (pa- 
Theorem 5.11 (Inversion theorem) Let ipa be the factorial umbra of a. It is 



g{t) = {i+tpr = j2{n)jpi-_ 



j=0 



with g.f. l + log[/(l + t)] by (31). 



with (5 the Bell umbra. 



Proof. By the Proposition 3.3 it results 



^a-X = Oi.x-P = a- 
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5.1 Factorial moments of the Poisson r.v.'s 

Being {'y.p.a).x = l-P-ifa) the factorial umbra of the umbral composition ^.(3.a is 
the umbral composition of 7 and the factorial umbra of a. Prom (13) the compound 
randomized Poisson r.v. Sn = ^1 + • — h Xjsi with N a Poisson r.v. with parameter 
the r.v. Y has factorial moments 

n 

gkBn,k[{ij)l, (m)2, • • • , {lj)n-k+l\ (34) 

k=l 

where (/x)i are the factorial moments of the r.v. X and are the moments of the 
r.v. Y . Now setting 7 = x.u in j.p.a, we have = x^. Then from (34) 

n 

x^Bn,k[{lj)i, {1^)2, • • • , {n)n-k+i] (35) 

k=l 

are the factorial moments of a compound Poisson r.v. with parameter x. Set a = u 
in x.p.a. We have {x.f3.u).x = x.P-x = x.u so that the factorial moments of x./3.a 
are equals to as well as for its probabilistic counterpart, the Poisson r.v. with 
parameter x. 

At the end set a = u in j.p.a. We have {'j.p.u).x = l-P-X = 7 so that the 
factorial moments of 7./? are equals to the moments of 7. Then a randomized Poisson 
r.v. with parameter a r.v. Y has factorial moments equal to the moments of the 
r.v. Y. 
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